Series for the Wave Function of a Radiating Dipole at the 

Earth's Surface 

By S. O. RICE 

In this paper three series expansions are derived for the wave 
function of a vertical dipole placed at the surface of a plane earth. 
Two convergent series and one asymptotic series are obtained. A 
remainder term for the latter series is given which enables one to set 
an upper limit to the amount of error obtained by stopping at any 
particular stage in the series. 

Introduction 

THE wave function above the earth of a vertical dipole placed at 
the surface of a plane earth is x 

*<r, .) - OV + M f . pg^gp , (i) 

where r and z are the horizontal and vertical distances from the dipole. 
k\ and fo are constants depending upon the electrical properties of the 
air and ground, respectively. 2 We shall be concerned with the value 
of this function at the surface of the earth. Setting z = gives us 
an integral for IIi(r, 0) which is the function of r to be investigated 
here. 

Although the electric and magnetic intensities are the properties of 
an electromagnetic field which have the greatest physical significance, 
writers on this subject often deal with the wave function because of 
its simpler form and because in many cases of practical interest it is 
nearly proportional to the electric intensity. However, the electro- 
magnetic field may be obtained from the wave function by differ- 
entiation. If the real parts of He~ iat and Ee~ iat represent the electric 
and magnetic intensities the field above the earth produced by the 
dipole is 

v dr 

ic> d'mjr, z) ic* 1 d I dUxir, z) 



oj drdz v u r dr\ dr 

1 A. Sommerfeld, Ann. der Physik, vol. 28, pp. 665-736, No. 4 (1909). 

2 The symbols used here are denned in a list at the end of the paper. 
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From these expressions it will be observed that if we obtain ex- 
pressions for IIi(r, 0) we shall be able to compute the field at the 
earth's surface except for the radial component E T , which is small 
compared to E,. 

Statement of Results 

The asymptotic expression for IIi(r, 0) is 



Mr, 0) = - 



1 



n = l 



(1 - r 2 )r . 

where Rm and R20 satisfy the inequalities 
(N + 1) I e^Vcsc"? 



e Z- — /.- Nw r -Kin — rKai , 



(irsr) n 



(9) 



\Rur\ < 



[r(ki - s) sin 0Y +1 



I-R20I < 



giftjr 



rko — rs\' 



6 = tt/2 — arg (ki — 5) being an angle slightly greater than tt/2. 
The convergent series for IIi(r, 0) are 3 



n »( r » °) = rr-^-r\rr\ eiklT £ C- isrr)»PL%r(.ki/s) 

\>- — T )T \ A |_ n _o 



sr 



- re**' £ ( £ PUfcWs) 

n=0 \*7" 



(14) 



and 



Mr, 0) = (^ [ £ ^ W. - nil ; 1/2 ; W ) 

(*jfe 8 r) n 



-r 2 E 



M! 



F(l, - w/2; l/2;s 2 /fei 2 ) 



(19) 



The quantities r and 5 are defined by t = fa/fa and 1/s 2 = 1/fci 2 + I/&2 2 , 
and the numbers on the right are the equation numbers in the text. 
W. H. Wise 4 has obtained series which are equivalent to those ap- 
pearing in (9) and (14). 

Procedure 

The results given here depend upon a transformation of the integral 
obtained by setting z = in equation (1). This integral can be 
expressed in the following way as has been shown by B. van der Pol : 5 



Hi(r, 0) = - ,-I-j | e — d{w> - i)-w, 
1 T Jk,i» 



(2) 



1 The Legendre functions are discussed by E. W. Hobson, "Th. of Spherical and 
Ellipsoidal Harmonics." Hypergeometric functions are discussed in Chap. XIV, 
"Modern Analysis," by Whittaker and Watson. 

* W. H. Wise, Proc. I.R.E., vol. 19, pp. 1684-1689, September 1931. 

6 Jahrbuch der drahtlosen Telegraphie Zeitschr. f. Hochfrequenz Techn., 37 (1931), p. 
152. 
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which becomes, after integration by parts, 



= + 



1 - T 2 






(3) 



The path of integration is the straight line in the complex w plane 
joining the points ki/s and fa/s. Arg (w — 1) and arg (w + 1) are 
taken to be zero at the point this contour crosses the real axis. The 
Argand diagram for a typical case is shown in Fig. 1. From the 
definitions of ki, kz, and 5 it follows that \s\ < | ki \ < \kt\, and 
= arg ki < arg s < arg ki < 7r/4. 



w = -h ISF 
= k|^ is A 
= k 2 / s is B 




Fig. 1 — Paths of integration in the w plane. 

Asymptotic Expansion 

To obtain an asymptotic expansion for IIi(r, 0) we deform the linear 

path joining A and B into the path ACDB as is shown in Fig. 1. 

The lines AC and BD are both inclined to the real axis at the angle 

arg (is*) where s* is the conjugate of s. This is the direction in which 

the exponential term e" rw decreases most rapidly since along it the 

variable part of the exponent is real and negative. 6 The section CD 

may be displaced to infinity where its contribution to the value of 

the integral becomes zero because of this exponential decrease. 

8 To show this for the line A C we set w = ki/s + is*u. As w goes from A to Cu 
is real and increases from zero. The exponent then becomes isrw = ik\r — |s| 2 rtt 
since ss* = \s\ 2 . 
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The integral IIi(r, 0) is then composed of two components consisting 
of the integrals along AC and DB, respectively, and we may write 

Mr, 0) = - (1 _ T r2)r [/(*0 - /OWJ (4) 

where 

J<\ 00 is* 
e i~*d(yfi - l)"l/2. (5) 

*/s 

We integrate (5) by parts N times and find 

[N ( — \rifjn -\totl* 

-'""iHri:^- 1 )" 1/2 

J roofs' ^iario fJN+l 
u pp.*' " «-** 



The derivatives may be expressed in terms of Legendre polynomials 
by means of the relation 

i-y4~ (™ 2 - I)" 1 '' = »K^ 2 - l)-" /2 - 1/2 Pn ( , = V 

K ' dw» K \Vw» - 1 / 

When the limits in the integrated portion are inserted and the definition 
of 5 used we see that 

where 

*- - - TBSF'JL, ^ L VP^l J (- 2 - D™ (7) 

An inequality for Rut may be obtained by using 3 

iiv+i(oi ^ i* + vr=ii*« 

which holds for all values of t in the t plane cut from — 1 to + 1, 
if arg V/ 2 — 1 = when * is real and greater than + 1. For then 
the absolute value of the Legendre polynomial in the integrand is 
seen to be less than \(w + l)/(«> - 1)| (JV+1)/2 when R(w) > 0, and R iN 
may be compared with an integral having \e iarw \ and powers of the 
factors | w + 1 1 and \w — 1| in the integrand. On the path AC we 

3 E. W. Hobson, loc. cit., p. 60. 
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have \w — 1| ^ \(k\/s — 1) sin 0| where 

9 = arg is* - arg ( -^ - 1 ) = ^ - arg {k\ — $) > -x • 



x . 2 « 8 v«i -, - 2 



Similarly we have | w + 1 1 ^ | fe/s + 1 1 • These inequalities enable 
us to deal with the integral of \e ilrw \ which may be integrated to 
show that 



\Rw\ < 



(N + 1) !e*Wcsc 6 
[r (Jfei - s) sin 0]" +1 



(8) 



By interchanging &i and fc 2 in (6) and (7) we obtain expressions for 
I(k%) and RiN. An inequality for \Rw\ is obtained from (8) by 
setting d = tt/2 and interchanging h and k 2 . By combining these 
expressions in accordance with equation (4) we obtain an asymptotic 
expansion for IIi(r, 0). 

In general, I(k 2 ) is negligible in comparison with I(ki) because k 2 
has a positive imaginary part which causes e ik * T to decrease rapidly. 
Since I{k 2 ) = R 20 , -&20 being the remainder after zero terms, we may 
obtain an inequality for I(k 2 ) by setting N = 0, 6 = v/2, and inter- 
changing h and k 2 in (8). Then from (4) we have the result 

1*. o) = - a^K?/-^ + *• - *] ■ » 

where Rm satisfies the inequality (8) and |i? 20 | < \e ik * r /(rk 2 - rs)\. 
Series for rii(r, 0) in Ascending Powers of r 



Put 






(10) 



and define K(k 2 ) as being obtained from (10) by interchanging ki 
and k 2 . By referring to equation (3) we see that / may be written 
in the form 



n i('- 0) = (1 j T2)r CJf(*i) " ^TOOl (n; 



We write 



L " ^2 £i(« — 1)! J t Vw 8 - 1 J' 

the infinite series being uniformly convergent. 



— 
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From Hobson's contour integral definition 3 of P n m (t) it can be 
shown that if R(m) < 1/2 

*=' ■« " V* r(l/2 - *») J, V^-l *"' 

where arg (w — 1) = <p, arg (w — — — t + ip, where <p is the angle 
measured counter-clockwise from the positive direction of the real axis 
to the line directed from w = 1 to w = t. Setting m = -f- 1/2 — n 
where n is a positive integer we obtain 

r (-W-_f)^ l dw = ( _ )n - 1(w _ 1} | Jk ( , 2 _ l)nl2 - m pl,2-n it y 

J x -yjw 2 — 1 \ 2 

Thus Equation (12) becomes 

r-r- (13) 

where in passing from the first to the second line we have set n = 
in 3 

***» " W+i) ( FTT ) "' ! "" , "(W 1/2 ; « + 1/2 ; if-'), 

and hav e summed the resulting series to show that P-ii 2 (ki/s) 
= V2/7TT. The function K(k 2 ) may be obtained from (13) by inter- 
changing k\ and k%. 

Combining (13) and (11), and using t = &1/&2 gives the convergent 
series for IIi(r, 0) given in the statement of results as equation (14). 

Another Power Series for I 

Here we obtain an expression for 7" somewhat similar to the one 
obtained in the previous section. The first step is to deform the 
contour joining the points A and B (w = kjs and w = k^s). The 
deformation is carried out in two steps shown in Figs. 2a and 2b, 
respectively. 

In Fig. 2(a) the contour joining A to B has been pulled around the 
point + 1 and looped over itself. The point II is destined to move 

3 E. W. Hobson, loc. cit., p. 188. 
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over to B and G is to move over to A. This deformation of the 
contour does not alter the value of the integral as long as we pay 
attention to the arguments of w - 1 and w + 1. In Fig. 2(6) the 
deformation is almost completed ; all that remains is for G to coincide 
with A and H to coincide with B. 





Fig. 2 — Deformation of contour in w plane. 

Using this deformation of the contour we may write equation (2) 
as follows: 



with the understanding that arg w — 1 and arg w + 1 have their 
principal values at the beginning of each integration. Upon referring 
to (2) we see that the middle integral is — IIi(r, 0) and hence 



where 



Ul{r ' 0) = 27(1 - r») [L( *° " L(fe) l 






(1+) e i>rw W( l w ^ « (^ f )» /'(1-f) W+l^fo, 

/fci/s k-i)«-£tt| i/s (^ - l)»/« 

and L(£ 2 ) is obtained from L{k\) by interchanging &! and k 2 . 



(15) 



(16) 



Let w°- - 1 = r 2 (l - /), or ot = Ml - ts 2 /k* 2 , then 

/.d+) w n+i dw __ &_ / M b r a+) [d - (s W)] n/2 

J*,. (^ 2 -D 3/2 2*xU/ J 



(1 ~t) 



(It 



(17) 



= 2|(^) n F(l,-n/2;l/2;5W) I 
where it is understood that at the initial point of the contour 
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arg (1 - /) = 0, arg (1 - (*W)0 = °- This ma y be verified b y 
expanding the numerator of the integrand and using 

r (1 +> v , . „ ■ v Vim + l)Y(v + 1) 

where w is a positive integer or zero, with v = — 3/2. 
Expression (16) now becomes 

Ukx) = 2%± ^- F(l, - n/2 ; 1/2 ; sW) (18) 

«1 n=0 » ! 

and the series converges for all finite values of r since the series inte- 
grated termwise in Equation (16) is uniformly convergent. 

We obtain the series for IIi(r, 0) given in statement of results as 
equation (19) by putting (18) and the corresponding expression for 
L(k 2 ) in equation (15). 

Notation 

The following symbols are used. C.G.S. electromagnetic units are 
used throughout the paper. 

c = velocity of light, 3 X 10 10 cm./sec. 

F(a, b;c;x) = The hypergeometric function 

ab a(a + l)ft(ft + 1) ■ ■ ... 

1+ T£* + 2!e(c+l) + - 

J (£r) = Bessel function of the first kind, zero order. 
ki = ufc. 

k 2 = Vcoj 2 + iliraa). The real and imaginary parts are 
positive. 

P n {t),PV.\y 2 n {t) = Legendre's polynomial, and associated Legendre's 
function of the first kind. 

Rin, Rn = Remainder terms in asymptotic series. 

r = horizontal distance of representative point from 

dipole. 

5 = fciWVfci 2 + W or 1/s* = 1/fci 2 + 1/M The real 
and imaginary parts of s are positive. 

s* = The complex conjugate of s. 
t = time in the introduction, otherwise a complex 
variable. 

w = complex variable. 
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2 = height of representative point above ground. 

e = dielectric constant of the ground in e.m.u. The 
dielectric constant of air in e.m.u. is 1/c 2 . If 
the dielectric constant in e.s.u. is e', then 
e = e'/c 2 . The dielectric constant of air in 
e.s.u. is 1. 

Hi(r, z) = Wave function for z & for a vertical unit dipole 
centered at the interface between air and 
ground. The wave function for a unit dipole 
wholly in air is obtained by multiplying the 
wave function given here by 2/(1 + t 2 ). By 
a unit dipole is meant the system obtained by 
letting the length / of a conductor approach 
zero while the current in the conductor ap- 
proaches infinity in such a way that II = unity, 
where the current equals the real part of Ie~ io " 
and does not vary with position along the 
conductor. 

IIi(r, 0) = Value of wave function at earth's surface. 

a = conductivity of the ground in e.m.u. If the con- 
ductivity is a' mhos per meter cube then 
a = lO" 11 *'. 

T = ki/k 2 . 

o) = angular velocity, radians/sec. 



